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Abstract 

We discuss the role of kaon-nucleon and nucleon-nucleon correlations in kaon 
condensation in dense matter. Correlations raise the threshold density for 
kaon condensation, possibly to densities higher than those encountered in 
stable neutron stars. 
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The possibility of kaon condensation in dense matter was suggested by Kaplan and Nelson 
[f|], and it has been discussed in many recent publications. IM- The basic idea is that the 
energy of a Kex, is lowered by interaction with nucleons, and in neutron star matter in 
beta equilibrium one expects negative kaons to be present if the energy to create one in the 
medium is less than the electron chemical potential, /i e , the energy required to add an extra 
electron. 

One complication in determining the low energy interactions of kaons with nucleons is 
that the interactions of kaons with protons in vacuo is dominated by the A(1405) resonance. 
Fortunately this does not affect the K~ interaction with neutrons, which is expected to 
dominate in neutron stars. In dense matter the effects of the A(1405) are likely unimportant, 
because the kinematic region of interest is far from the resonance. To describe the KN 
interactions Brown, Lee, Rho, and Thorsson ||, herafter denoted by BLRT, employed an 
effective Lagrangian based on chiral perturbation theory. To circumvent the difficulties with 
the resonance, parameters of the Lagrangian were fitted to K + N, (N = n or p), data, which 
are unaffected by the resonance. The Lagrangian was treated in mean field theory, and 
the density for condensation was found to be ~ 3 — 4 p 0 , where p 0 = 0.16fm -3 is nuclear 
matter density. This is to be compared with the central density of about 4 p 0 for a neutron 
star of mass 1.4M 0 according to the estimates of Wiringa, Fiks and Fabrocini § using 
realistic models of nuclear forces. The conclusion was that kaon condensation could affect 
the structure, maximum mass and cooling rates of neutron stars significantly. 

In this Letter, we first show that an important contribution to the kaon energy in the 
BLRT treatment is the first term in an expansion in powers of the particle density, with a 
strength proportional to kaon-nucleon scattering lengths. We then argue that kaon-nucleon 
correlations in the medium will be reduced compared with those for a kaon and nucleon in 
vacuo , and that consequently the attraction experienced by a kaon in matter will be smaller 
that earlier estimates indicated. 

BLRT use the Lagrangian of Jenkins and Manohar [Jl] containing terms of the two lowest 
orders in the chiral expansion, and the nucleons are treated in mean field theory, that is 
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their densities are assumed to be uniform, and equal to p n and p p . The spatially uniform 
state, which is the one of lowest energy in this model, is given by 


K ± = v K e ±i€Kt 


( 1 ) 


and 


C = (1 - a)k + k~ - (m 2 K - 5 2 )K + K~ 
- l -V{K~k + - K~K + ), 


( 2 ) 


with V = (p n + 2p p )/(2/ 2 ), 5 2 = [T, KN (p n + p p ) - C(p n - p p )]/f 2 , and a = ~[D(p n + p p ) - 
D\p n — p p )]/f 2 - The V and 5 2 terms describe interactions of kaons with effective vector 
and scalar fields, while a takes into account energy dependence of the interaction. The pion 
decay constant, /(~ 93 MeV), and the constant C(~ 33.5 MeV) are well determined, but 
the values of T,kn, D, and D' are not well known. The experimental K + N scattering lengths 
Toll are used to determine D for chosen values of from Eq.(27) of BLRT. Eq.(28) of 


BLRT gives D' = 0.092 /rriK- A large T,^n implies a large energy dependence for the K~N 
interaction, and the value T,kn — 278 MeV was considered to be realistic by BLRT. It is 
presently considered likely that Y*kn lies around 400 MeV ||. The corresponding values of 
D are —0.23/mA' for E/v/v = 278 MeV, and — 0.48/mx for Y^kn = 4R0 MeV. 

The K~ energy may easily be calculated from Eq. (2), and is 
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as used by BLRT. Observing that a, 5 2 and V are linear in the nucleon densities, p n, we 
expand ex in powers of p N . To second order in the densities one finds 
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The first term is just the kaon rest mass. The second term may be expressed simply by 
making use of the fact that the K~N scattering lengths [jllj are given by 
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Q>k~p — —79- ( m K + ^kn + C + (D + .D )m 2 K ), ( 6 ) 

Airj-rriK 

where m R = miymK /(tun + mjf) is the reduced mass. This term is thus the Lenz potential, 

,Lenz (Pri®A' - n T PpOjR~p) 1 (7) 

m R 

the well-known result for the energy shift to first order in the density ( see also Ref. j7|). 
Thus for given values of the scattering lengths, the Le n z potential is independent of the 
choice of T, KN . However, terms of order p 2 N and higher depend on Tj KN . When T, KN is 
small, 5 2 and a are small (see Table 1) and the higher order terms dominated by V 2 /(8i71k) 
are repulsive. In contrast, when T*kn is large, 5 2 is large, 1 — a can become zero at relatively 
small densities, ~ lOpcb and the higher order terms are attractive. 

The dependence of the kaon energy on nucleon density is illustrated in Fig.l for a proton 
fraction of 0.1. The curves labelled by T, KN values are obtained from Eq.(|). For comparison, 
the Lenz approximation, which contains only the first two terms in Eq.([|), is also shown. The 
Lenz term gives a major contribution, while other terms representing the energy dependence 
of the interaction and relativistic effects may be significant. 

The Lenz term is exact only at low density, and there are corrections to it when it is no 
longer a good approximation to assume that the K~N relative wave function, i^k-Ni is th e 
same as for a pair of particles in vacuo. For s-wave scattering, if} at low energy is given by 

ipK-N(r) = 1 - a K - N /r for r > R, (8) 

where R is the range of the K~N interaction. The corresponding unperturbed wave function 
is simply unity. 

To illustrate the effects, we represent the K~N interaction by potentials, corresponding 
to the two possible isospin values, / = 0,1. The scattering lengths are given by 

O-K-N — f ( ^ 3r '' /? A-Ar(' r )V ; A“w( r ), (9) 
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where Vr-r ( r ) represents the K~N potential. Since 'i/’r-r/t <R) > 1 — or-r/R when 
Vr-n(t) does not have a repulsive core, 2ira K - N /mR is larger than the volume integral of 
Vr-n{t) by a factor larger than 1 — a K - N /R. The potential Vr-r(t) is believed to be short 
ranged, because the vector and scalar mesons that contribute to it have masses larger than 
the inverse of hadron radii. Hence Vr-r has a range R < 1 fm, and the factor 1 — cI r-r/R 
is larger than 1.41 and 1.85 for N — n,p with the scattering lengths dR- n = —0.41 fm and 
a,R- p = —0.85 fm given by Eqs.(5) and (6), and the parameters taken from BLRT. These 
factors represent the correlation between low-energy K~ and nucleons in vacuo , and when 
it is significantly different from unity, the Lenz term is valid only at densities low enough 
that the radius, r 0 , of a sphere with a volume equal to the average volume per particle, 
1/p = 47rrQ/3 is larger than both R and |a|. 

In order to make quantitative estimates, we fitted square well and Yukawa potential to 
the scattering lengths. For a square well of radius 0.7 fm, the depths are 126.2 MeV for 
neutrons and 181 MeV for protons, while for R=1.0 fm, the corresponding values are 49.5 
and 75.6 MeV. For a Yukawa potential V 0 exp(—m x r)/(m x r), one finds V 0 = —795 MeV for 
neutrons and -1075 MeV for protons if the mass, m x , of the exchanged particle is taken to be 
that of the p meson, 770 MeV. In neutron stars, it is the K~n interaction that is the more 
important, since the proton fraction of matter is small, and in Fig. 2 we show ipR- n for the 
three potentials. These exhibit the enhancement of the wave function for small separations. 

At higher densities, the kaon energy in matter should be calculated either by summation 
of diagrams, as in approaches that evaluate the K~N interaction G matrix and higher-order 
correlations, or from variational wave functions of the type 

T = \n i=liA f K - N (r K - ~ ri)]Tiv(fi, (10) 

where J'r-r represents the K~— nucleon correlation and Tjyis the ground state of a system 
with A nucleons. At the large densities of interest in the cores of neutron stars, ro < 1 fm. 
When r'o < R, , the kaon is weakly correlated with nucleons, and .I’k-r is unity to a good 
approximation. In this limit the Lenz potential is replaced by the Hartree potential, and 
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the K energy is given by 


{-Hartree = % + J [p n V K - n (r) + p p V K - p {r))d 3 r. 


( 11 ) 


The Hartree potential is less attractive than the Lenz one by a factor in excess of 1 — 
O’K-n!R- In Fig- 1 we show estimates of the Hartree contribution to the kaon energy for the 
potentials we fitted to the scattering lengths, as described above. The Hartree calculation 
does not include effects of energy dependence of the interaction or of relativity, but since the 
Hartree energy is closer to mx than is the Lenz term, we expect these effects to be smaller 
than estimated from Eq. (|4| ). 

To investigate the validity of the Hartree approximation we now estimate second order 
terms in the K~N interaction. For simplicity, we neglect the effects of the small number of 
protons. In this case, the main coupling of the kaon is to neutron density fluctuations, since 
the kaon is spinless, and consequently does not couple to spin fluctuations. The contribution 
to the energy of a kaon at rest is then given by 




( 12 ) 


Here < l\pq\0 > is the matrix element of the operator for the Fourier transform of the 
neutron density between the ground state, 0, and an excited state, l, and E 0 and Ei are 
the energies of the states. The energy of a kaon with momentum p is denoted by £k(p)- H 
one neglects the recoil energy of the kaon, one finds an upper bound for the energy shift in 
second order: 


where 


- Ae K < \ Y. I v 0?) \ 2 x(q, u = 0) 

Q 

= -/ d 3 rd 3 r'V(r)x(r — r')V(r') 


X(q, 0 ) = £ 2 

i 


< l\ Afl° > 
Ei — E 0 
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(13) 


(14) 
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is the static density-density response function for neutrons, and y (r) is its Fourier transform. 
The recoil energy of the kaon is comparable in magnitude to the excitation energy of the 
neutron liquid for the wavevectors of importance, so the actual second order contribution 
to the energy is expected to be considerably smaller than the bound. The long-wavelength 
static density-density response function for neutron matter may be estimated by using the 
thermodynamic identity x(q — 1 > 0, u; = 0) = dp n /dp , n , together with many-body calcula¬ 
tions of the neutron chemical potential as a function of density. The dependence of x on 
wavenumber is not well known, and to obtain an estimate of the second-order term we shall 
replace x(q) by its long-wavelength value. At a density of 4p 0 , the estimates of the second 
order term obtained by using values of dp n /dp n from Ref. [12[ are 7 MeV (square well, R = 1 
fin), 16 MeV (square well, R = 0.7 fm), and 46 MeV (Yukawa), and for higher densities 
the estimates are smaller. In all cases, these energies are considerably less than the Hartree 
potential. 

The second order contribution to the energy is reduced greatly due to the strong repulsive 
interactions between nucleons. Had we neglected neutron-neutron interactions, the estimates 
of the second order term, which would then reflect Pauli blocking of intermediate states 
in the repeated scattering of a kaon by a neutron [13], would have been much greater. 
For a free Fermi gas, p n dp n /dp n = 2Ap/3, where Ep ps 58(p n /po) 2 ^ 3 MeV is the Fermi 
energy. At a density of 4p 0 , p n dp n /dp n estimated from the calculations of Friedman and 
Pandharipande [12[ is about 5 times larger than the free gas value. The fact that stable 
neutron stars with a mass of 1.4 M e exist is a clear indication that the equation of state of 
neutron star matter is considerably stiffer than that for a free Fermi gas. While there are 
still uncertainties about the properties of matter at high densities, this increased stiffness is 
a rather general property shared by many models, as one can see from, e.g., Ref. [|[]. As a 
first step towards a more complete many-body treatment, we are calculating the kaon energy 
using the lowest order constrained variational (LOCV) method [[L4j, and the results will be 
reported elsewhere. 

To explore consequences of our calculations for kaon condensation, we also show in Fig.l 
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the electron chemical potential, /i e , calculated by Wiringa et al. 0. Our Hartree estimates of 
kaon energies always lie above /i e for essentially all stable neutron star models. However, for a 
number of reasons we are prevented from definitively ruling out kaon condensation in neutron 
stars. First, we have used simple approximate forms for the kaon-nucleon interactions, and 
estimates should be made with better forms of the interaction, which should also take 
into account the variation of the coupling constant, /, with nuclear density ||. Second, 
corrections to the Hartree result need to be investigated in detail. Third, as one can see 
from the figure, estimates of ji e at high densities are subject to considerable uncertainty: 
models of nucleon interactions that fit laboratory data equally well lead to very different 
values of /i e in neutron star matter. All these points need to be looked into. 

The main conclusion of our calculations is that, in a dense medium such as that in the 
interiors of neutron stars, correlations are very different from those for a dilute gas. This 
has the effect of reducing the attraction experienced by a kaon in dense matter. While a 
low-order chiral expansion may be useful for describing low-energy scattering of kaons by 
nucleons, this information alone can predict the kaon energy only at low densities. When 
the K~N scattering lengths are as large as indicated by BLRT, the Lenz approximation 
overestimates the attraction felt by a K~ in neutron star matter at the density of neutron 
star cores by roughly a factor of 2. 
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Nuclear Theory at the University of Seattle for hospitality during the workshop there on 
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FIGURES 


FIG. 1. Energy of a single K~ in matter as a function of density, and electron chemical po¬ 
tentials taken from the calculations of Ref. [6]. The calculations of kaon energies are for a proton 
fraction of 0.1. Curves labelled by values of T^kn, hr MeV, are evaluated from Eq.(3). Hartree 
potentials, evaluated from Eq.(ll), are shown for square wells of radii R = 1 fm and R = 0.7 
frn, and for a Yukawa potential. A and U refer to the AV14 and UV14 two-body potentials, and 
VII and TNI to the UVII and TNI many-body potentials. The open circles indicate the central 
densities of neutron stars with a mass of 1.4M 0 , and filled circles indicate the central density for 
the most massive stable neutron star for the given equation of state. 

FIG. 2. Wave function for the relative motion of a K~ and a neutron at zero energy for square 
well potentials (solid lines) and a Yukawa potential ( long-dashed line). The short-dashed line 
shows the asymptotic limit of the wave functions. 
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Fig.l, Pandharipande et al., Kaon Energies in Dense Matter 





